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ABSTRACT 

We demonstrate that pulsar timing measurements may potentially be able to detect 
the presence of dark matter substructure within our own galaxy. As dark matter 
substructure transits near the line-of-sight between a pulsar and an observer, the 
change in the gravitational field will result in a delay of the light-travel-time of photons. 
We calculate the effect of this delay due to transiting dark matter substructure and 
find that the effect on pulsar timing ought to be observable over decadal timescales for 
a wide range of substructure masses and density profiles. We find that transiting dark 
matter substructure with masses above 10~^ Mq ought to be detectable at present by 
these means. With improved measurements, this method may be able to distinguish 
between baryonic, thermal non-baryonic, and non-thermal non-baryonic types of dark 
matter. Additionally, information about structure formation on small scales and the 
density profiles of galactic dark matter substructure can be extracted via this method. 

Key words: gravitation, pulsars: general. Galaxy: halo, dark matter, large scale 
structure of Universe 



1 INTRODUCTION 

One of the most fundamental goals of cosmology is to deter- 
mine the composition of the universe. Contrary to our ex- 
perience in the solar system, where the emitted light traces 
the mass distribution, the vast majority of matter in the 
universe appears to be non-luminous. This puzzle is known 
as the dark matter problem. One of the greatest goals in 
modern cosmology is to determine the nature and proper- 
ties of this dark matter. Although many solutions have been 
proposed explaining the presence of dark matter, its nature 
remains mysterious, as all but the most indirect methods of 
detection have been unsuccessful thus far. 

The astrophysical evidence for the existe nce of dark 
matter is overwhelming, as recently reviewed in lSiegel et al.l 
l|2006l ). While the amount of matter in luminous sources is 
only n* ~ 0.005 (|Fukugita fc Peebles! [20oi : IShankar et al.1 
I2OO6I ). precis ion measurements of the cosmic microwave 
background jSpergel et al.l I2OO6I ) and type la supernovae 
less et al]l2006l ) indicate that the total amount of mat- 
ter is Q.m — 0.26, or that 98 per cent of the matter in 
the universe is non-luminous. While a significant compo- 
nent of this non-luminous matter can be baryonic, nu- 
cleosvnthesis llFiel ds fc Olivel lioO^ . Silk damping effects 
iMiller. Nichol. fc Chen .2003 ). an d searches for MA ssive 
Compact Halo Objects (MACHOs) (|Alcock et al.ll 19981 ) con- 



strain the amount of baryonic matter to be 57;, ~ 0.04. This 
means that, although a portion of the dark matter is bary- 
onic (~ 13 percent), the vast majority of matter in the uni- 
verse is both non-luminous and non-baryonic (rioM). This 
non-baryonic dark matter must also be cold; i.e., it must 
have become non-relativistic at early times, as evidenced 
most clearly by the large amount of power on small scales 
at early times in the Lyman-a forest ( Viel et a l. 2005). This 
constraint eliminates the only standard-model candidate for 
non-baryonic dark matter, the neutrino. There are then two 
generic classes of candidates for cold, non-baryonic dark 
matter: the thermal relic and the non-thermal relic. A ther- 
mal relic is a particle produced by thermally efficient inter- 
actions in the background plasma of the early universe. Its 
abundance then freezes out as the temperature drops, but it 
remains in kineti c equilibrium with the background plasma 
until later times ('Berts chingeij fiooil ) . Examples of thermal 
relic dark matter candidates include the lightest supersym- 
metric particle, lightest Kaluza-Klein particle, and other 
particles generically classified as Weakly Interacting Mas- 
sive Particles (WIMPs). On the other hand, non-thermal 
relics can be produced at phase transitions by vacuum mis- 
ahgnment mechanisms, and are born non-relativistic. Exam- 

! iles of non-thermal dark matter candidates include a xions 
Peccei fc Quinnl ll977': 'Weinbere"l978': 'Wilczek^'l978^ and 
massive gravitons |Dubovskv, Tinyakov, fc Tkachcv 2005). 
Finally, it is possible that both thermal and non-thermal 
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types of dark matter exist, and that JIdm is made up of a 
composite of these two classes. 

In a recent and exciting development, non-baryonic 
dark matter has been detected indirec tly via its gravita - 
tional influence in the Bullet Cluster l|Clowe et al.l l200^ . 
Since baryons, thermal relic dark matter, and non-thermal 
dark matter all decouple from the primordial plasma at 
different times and possess different interactions, it is ex- 
pected that the small-scale structure that forms composed 
of these types of matter will be quite different from one 
another. Baryons will remain in either diffuse gas clouds or 
collapse to form MACHOs, thermal dark matter will have its 
small scale structure on scales below 10"* to 10~® Mq sup- 
pressed by Silk damping, whereas non-thermally produced 
dark matter will have no such cutoff in its spectrum. Ther- 
mal dark matter will produce, on small scales, diffuse WIMP 
microhalos after gravit ational collapse down to masses of 
about one Earth mass (iGreen. Hofmann. fc Sc hwarz 2005; 
Loeb fc Zaldarriagall2005l:lDiemand. Moore, fc Stadel.2005i: 



Bertschin gcr 2006), whereas non-thermal dark matter will 



produce collapsed structure down to much smaller mass 
scales, without a cutoff in the spectrum, which will produce 
Non-thermal Axionic Collapsed HalOS (NACHOs). Under 
certain circumstances l|Sikiviel2006l : IZurek. Hogan. fc QuiiinI 
l2006l ). the non-thermal dark matter will produce halos 
with substantial enhancements in density and abundance 
at ~ 10"^'^ Mq, known as axion miniclusters. In either case, 
these dark matter structures on the smallest scales may sur- 
vive intact to the present day as substructure within larger 
collap sed structure such as our own galaxy (jCoerdt et al.l 
[2006). While large-scale dark matter structures like our 
galaxy's halo are expected to be smooth to a simple approxi- 
mation, realistic models include clumps, cusps, and possibly 
caustics. 

This paper demonstrates that pulsar timing measure- 
ments can be used to probe the dark matter substructure 
in our own galaxy which transits across the line-of-sight be- 
tween a pulsar and an observer. As the dark matter sub- 
structure transits, the variation in the light-travel-time of 
the pulses is of sufficient magnitude that the more stable 
millisecond pulsars should be able to not only detect the 
substructure, but to discriminate between the different sig- 
natures of MACHOs, NACHOs, and WIMPs. The remain- 
der of this paper is layed out as follows. Section [2] presents a 
derivation of the delay in the light-travel-time due to a grav- 
itational source intervening near or across the line-of-sight 
between an emitter and an observer, providing a number of 
examples. Section [3] applies this idea to the physical case of 
using pulsar timing measurements to detect the intervening 
dark matter substructure between Earth and the pulsar of 
interest. Section |3] discusses the possibility of differentiating 
between different types of dark matter and constraining its 
properties through precision timing measurements. Finally, 
Section [S] summarizes the major successes and drawbacks of 
this method, detailing the challenges of using pulsar timing 
to search for dark matter in our own galaxy. 



2 GRAVITATIONAL TIME DELAY 

One consequence of general relativity is the delay in the 
light-travel-time (l.t.t.) due to intervening gravitational 




Figure 1. The generic configuration for a fixed emitter and ob- 
server, with a gravitational source of radius rmax transiting across 
or close to the line-of-sight. The emitter-observer distance is taken 
to be Z, while the gravitational source has a distance to the ob- 
server X and a velocity v, both of which have components per- 
pendicular (xx, I'x) and parallel (xy, U||) to the line-of-sight. The 
minimum distance between the gravitational source and the line- 
of-sight as the source transits is defined as the impact parameter, 
b. 



sources between an emitter and an observer (|Shapirci|ll964l ). 
This effect, known (for a point source) as the Shapiro time 
delay, is o ne of the classic sol ar system tests of Einstein's 
relativity (|Baez fc Bunnll2005^ . This effect can be applied 
to extended sources and larger scales as well. If any col- 
lapsed substructure within our galaxy transits across the 
line-of-sight (l.o.s.) between an emitter and an observer, then 
the l.t.t. will be altered by the gravitational presence of the 
galactic substructure. The remainder of this Section is de- 
voted to a derivation of the magnitude of the delay in the 
l.t.t. due to transiting collapsed substructure. 

Consider the configuration shown in Figure I: an emit- 
ting source is a fixed distance I from an observer, with a 
transiting gravitational source near the l.o.s., with a min- 
imum distance to the l.o.s. (i.e. impact parameter) b. The 
gravitational source has a velocity v{t) and position x{t), 
where the velocity, assumed to be constant, has components 
parallel (u||) and perpendicular (vx) to the l.o.s., while the 
position deviates from its initial value only in the direction 
given by v. The two cases of physical importance to consider 
are when the radial size of the gravitational source (r-max) is 
either greater or smaller than the impact parameter, h. 

To determine the time delay induced by the presence of 
this gravitational source, we integrate a null geodesic from 
the emitter to the observer. We consider the weak-field met- 
ric for a single gravitational source. 



ds^ = -(I + 2^)^dt' + {1- 2 J)rfx' + X^d^\ 



(1) 



where (j> is the gravitational potential induced by the pres- 
ence of a single gravitational source. In the presence of mul- 
tiple sources, we note that this procedure can be used to 
calculate the time delay for each source individually and 
then the effects can be summed, since the fields are weak. 

We choose a radial null geodesic along the l.o.s. from 
the emitting source to the observer, working in coordinate 
time, and find that the total l.t.t. is 



t = to + M = 



dt = 



i(I-2^)dx. (2) 
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The first term on the right-hand-side of equation ((2)| inte- 
grates to give a constant (l/c) which we identify as to, or the 
i.t.t. in the absence of any gravitational sources. The second 
term is indicative of the time delay (At) due to the presence 
of a gravitational source. 

Since the gravitational fields of interest are weak 
(i^/c^ ^ 1), it is an excellent approximation to treat (f> as a 
Newtonian potential when outside the gravitational source, 
(j) = —GM/r. When one is inside a spherical gravitational 
source, one must be a little more careful and perform the 
following integral: 



,3 , Gp{r') 
a r 



(3) 



which reduces to (j) = ~GM/r at r = rmax- In general, the 
effect of the time delay can then be quantified by the integral 
from the emitter to the observer 



At ■ 



dx' , 



(4) 



where the r from equation ([3} is the distance between each 
given position along the l.o.s. {x') and the center-of-mass 
of the gravitational source. For points lying outside of the 
source, the total mass (Mtot) of the source contributes to 
(f), whereas the gravitational potential can either increase or 
decrease as one penetrates the source, dependent upon the 
density profile. 

For an extended source of radius rmax, the enclosed 
mass within the source is given by 



At ~ 2.95 X 10" 



In 



/42:||(/ - x{) 



(8) 



where is one Earth mass. As the position of the grav- 
itational source, x, changes with time, so too will the time 
delay. At. We note that the time delay of equations ([Hl-IHI 
depends very weakly on changes in I and x\\, but is more 
sensitive to changes in x±, since x± <C , Z — a;|| }. We also 
note that these results indicate that the time delay scales 
linearly with mass, but only logarithmically with position. 

As the gravitational sources are expected to be extended 
objects if they are non-baryonic, we also consider the case 
that the l.o.s. intersects the gravitational source, so that 
b < rmax. The general equation will still be given by equa- 
tion Q, although the integration will be significantly more 
difficult here. Generically, the integral must be broken up 
into three parts: from the observer to the near edge of the 
source, from the near edge through the source to the far 
edge, and from the edge of the source to the emitter. The 
integral of equation Q then becomes 



At 



+ 



-V('-max)2-(xi)2 



||+V('-max)2-(xi)2 



C-^]dx' 



-V('-ma,x)2-{l_L)2 



dx' 



C-^]dx' 
■ \ c^r 



(9) 



4Tvr p{r)dr, 



(5) 



but merely substituting Menc for Mtot is insufficient; one 
must recalculate <j) at each point and then solve equa- 
tion ((4|. The total mass of the gravitational source (Mtot) 
is given by equation ([5]) at the point where rmax = 
[x' — a;||)2 + {xa_Y . In general, therefore, the Shapiro time 
delay will be given by the combination of equation Q with 
the appropriate value for at all points. 

We first consider the case of impact parameters, h, that 
are larger than the maximal radius (rmax) of the gravita- 
tional source in question. This approximation corresponds to 
the physical cases of baryonic substructures (stars and MA- 
CHOs) and for non-baryonic substructure where the object 
is sufficiently far away, even at its point of closest approach, 
so that its radius never intersects the l.o.s. between the emit- 
ter and observer. For this case {h > rmax), the enclosed mass 
is always Mtot, and equation (|4]) is exactly integrable, yield- 
ing a time delay of 



_ 2GMtot 1^ r - ^11 + - ^ll)" + (^^) ] .g) 

In the region of physical interest, where the approximation 
{I, x\\,l — x\\} 3> {&, x±} is valid, a simpler expression of 



_ 2GMtot / 4a||(/-x||) 

{XA_Y 



(7) 



can be obtained. Numerically, this yields a time delay de- 
pendent only on the parameters I, x, and Mtot, such that 



where (j>{^') is given by equation (|3]), and is valid for any 
generic density profile p(r) for the gravitational source. An 
oftentimes more useful form for equation ([SJ is the partially 
solved equation 



At 



2GMt, 



■In 



4a:||(Z — a::||)(rmax — (rmax)^ — {x±)^) 

(a;x)^(rmax -I- \/ (nnax)^ — (x l_y^) 



II 



-md^r. 



V('-max)2-{a;_L)2 



(10) 



which contains the approximations x± ^ {x\\,l — x^^}, but 
still admits a generic density profile p(r) (implicit in <t>{r)). 
For all physical cases (i.e., p(r) is everywhere non-negative), 
this leads to a smaller time delay compared to the case 
of a point mass for sufficiently small values of x±. Fig- 
ure [2] illustrates At as a function of x± for the case of 
a point source, the case of a constant density halo with 
p(r) = 3 Mtot/47rrmay, as well as the case of a halo with 
an NFW (|Navarro. Frenk. fc Wh"it3l 19961 . [l997l ) profile with 
a cutoff at rmax, obeying the formula 



p{r) = 



Mt< 



47rr(ro + r)2 [ln(l + ^) 



(11) 



ro + ri. 



where ro is the turnover radius between an outer halo density 
falling off as and and inner halo density which falls as 
1/r. As can be seen in Figure [5J smaller impact parameters 
are more sensitive to the form of a gravitational source's 
density profile. 



© 2007 RAS, MNRAS 000, 000-000 



4 E. R. Siegel et al. 




0.0005 0.001 0.0015 0,002 

(pc) 



Figure 2. At vs. x±_ with I = 21{pc, = Ikpc and Mtot = 
10'* Mg) for tiiree sources: a point source (solid line), a constant 
density halo (dotted line), and an NFW profile (dashed line). 
Additionally, the physical size of the gravitational source for ex- 
tended objects, Tmax) is chosen to be fO~^ pc. The NFW profile 
illustrated has a turnover radius ro ~ 10~*'^ pc. x±_ is plotted in 
units of parsecs, At in units of seconds. This result scales linearly 
with mass, as indicated in equation (|8]l, so that a microhalo 100 
times as massive would have a change in time delay 100 times as 
large. 



3 DARK MATTER DETECTION 

The analysis in Section [2] has demonstrated that a gravita- 
tional source lying near the l.o.s. between an emission source 
and an observer will cause a delay in the l.t.t. of the emit- 
ted signal. We now specifically turn to the case of our own 
galaxy for a physical application. It is not possible to mea- 
sure the absolute l.t.t. of any astrophysical object directly, 
but changes in the l.t.t. can be observed over time. To mea- 
sure a change in the l.t.t., we require a stable, predictable 
emitter whose emission properties remain consistent over 
long timescales. The most stable, consistent emitters in the 
known universe are the millisecond pulsars, many of which 
have remained stable without glitches, starquakes, or flux 
unpinning over timescales exceeding thirty years. Pulsars 
have been used in many instances in t he past to glea n in- 
formation about the physical universe (|Detweileijll979l ). By 
treating millisecond pulsars as emission sources, the l.t.t. of 
each pulse as it heads towards Earth will be quite sensitive to 
any changes in the matter distribution near the l.o.s. Thus, 
sufficiently accurate measurements of the l.t.t. will allow us 
to search for all massive transiting objects near the l.o.s., 
including MACHOs, WIMPs, and NACHOs. 

Pulsars are the ideal candidates to serve as emission 
sources, as they are the most perfect natural clocks in the 
known universe. Pulsars exist with a wide variety of ro- 
tational periods and also in the emission periods of their 
pulses. The pulsars with the highest rotational frequencies, 
and hence the shortest pulse-to-pulse periods, also happen 
to be the most stable types of pulsars, with a period (Tp) 
of 0(1 ms). These millisecond pulsars have typical residu- 
als of 0(1 fis), meaning that if one measures enough pulses 
accurately (typically over ~ 1 hr timescales), the deviation 



of any pulse-to-pulse measurement from the average value 
is 0(1 ps). Remarkably, these residual uncertainties are not 
cumulative, so that the arrival-time uncertainty between the 
first pulse and the n-th pulse remains 0(1 /is), and does not 
increase as the number of pulses increases. Because of these 
properties of millisecond pulsars, the ones which are com- 
pletely stable (no glitches, starquakes, or flux unpinning) 
have their periods known to an accuracy of ~ 10~^^ s or 
better, as an 0(1 /is) uncertainty over a span of 0(10^^) 
pulses (the flrst millisecond pulsar was discovered in 1982) 
yields a measured period accurate to ~ 10~^* s. 

In order to detect the time delay induced by transit- 
ing substructure, the effect must, in some sense, be large 
when compared with the intrinsic timing uncertainties. The 
uncertainty in the pulse arrival time is primarily due to en- 
vironmental uncertainties, owing to the difficulty in mod- 
elling both our own and the pulsar's local environments. 
In principle, the simplest method for detecting substructure 
gravitationally would be to determine the shift in the ex- 
pected pulse arrival time (ATp) for adjacent pulses due to 
changes in the gravitational potential near the l.o.s. This 
can be computed by examining the shift in At over the 
time of one pulse, Tp, as the intervening matter changes 
position. Assuming a constant velocity for the intervening 
gravitational source, ATp is calculable in general from equa- 
tion (|10p by taking the difference in At between an initial 
source position Xi = (b,x\\,x^^ and a final source position 
Xf — {b, -t-uiiTp, Xz + v±Tp). So long as ATp is greater than 
the known period accuracy (~ 10~^^s), it may be possible 
to extract the signal that would be caused by transiting 
substructure using statistical techniques, despite the large 
pulse-to-pulse residuals. ATp as a function of time is shown 
in Figure [3] for the three transiting cases of a point mass, 
a constant density extended object, and an extended mi- 
crohalo with an NFW profile. Unfortunately, this method is 
not presently feasible in practice for the detection of masses 
significantly below 1 Mq , since pulse arrival times have not 
been measured to accuracies even approaching 10~^^s. 

A more complicated, but perhaps more feasible method 
would be to integrate the observed and expected time delays 
over many pulses. For this method, we are not constrained 
by the timescale Tp, as we are free to measure ATp over 
any time interval that is a multiple of Tp. This allows us to 
enhance the magnitude of the signal from the pulse-to-pulse 
values of Figure [3] up to the difference in the time delay 
between two given positions. The method of computing the 
change in the l.t.t. is the same as above, except Tp in the 
final position is to be replaced by nTp, where n is any pos- 
itive integer. While the noise residual does not decrease for 
this case as more pulses are observed, it does not increase 
either, as the 0(1 /is) uncertainty is a non-cumulative error. 
Therefore, as long as the change in the time delay over a 
reasonably measurable timescale is comparable to or greater 
than the residual uncertainties, this signal will be detectable 
in practice. Calculations for the integrated effect of the time 
delay over many years for a transiting point mass, extended 
constant density halo, and an extended halo with an NFW 
profile are illustrated in Figure 4. Quantitatively, for a resid- 
ual of 0(1 us), a transiting mass of Mtot ~ 10* with an 
impact parameter of & ~ 10"* pc will be at the limit of de- 
tectability. Detection is made more difficult by the fact that 
modern pulsar analysis techniques require the removal of lin- 
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Figure 3. ATp vs. time for a set of parameters corresponding to 
tile variabies deiined in Figure |2] I = 2ltpc, xy = Ikpc, Mtot = 
10* M®, b = 10"* pc, vj_ = 200km s-\ and Tp = IQ-'-^s. ATp 
is illustrated in seconds and time in years. The solid curve is for 
a point mass, the dotted curve is for a microhalo of constant 
density with size T'max = 10 pc, and the dashed curve is for an 
extended microhalo with an NFW profile of size T'lxiax = 10 pc 
and turnover radius ro = 10~*'^ pc. 
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Figure 4. At vs. time for a set of parameters corresponding to 
the variables defined in Figure |2] I = 2kpc, = Ikpc, Aftot = 
10'* M®, b = 10-* pc, vj_ = 200km s"!, and Tp = IQ-'^s. ATp 
is illustrated in seconds and time in years. The solid curve is for 
a point mass, the dotted curve is for a microhalo of constant 
density with size T'max = 10 pc, and the dashed curve is for an 
extended microhalo with an NFW profile of size T'jnax = 10 ^ pc 
and turnover radius ro = lO"*'^ pc. The chosen set of parameters 
is very optimistic, but demonstrates the potential power of this 
method. 

ear and quadratic terms from millisecond pulsar data before 
any further analysis is done. The majority of close dark mat- 
ter encounters (excepting actual observed transits) will have 
their signatures absorbed in this subtraction. We suggest, 
therefore, that it may be beneficial to develop non-standard 
analysis techniques to explicitly search for this effect. 



As can be seen from Figures [3] and 4, the time delay 
signal is very sensitive to the density profile of the tran- 
siting source, and can be used to determine the densities 
of sufficiently massive transiting dark matter, if their im- 
pact parameters are small enough. Halos with steeper den- 
sity profiles and smaller cores are more easily detectable, 
whereas the more diffuse ones will have a more significant 
departure from the point-source template. Of course, at suf- 
ficently large impact parameters, all source behave as point 
sources. Many of the data analysis techniques being pio- 
neered by groups searchin g for a cosmic gra vitational wave 
background with pulsars l|Jenet et al.ll200q ) should be ap- 
plicable to identifying the time delay signature of transiting 
dark matter substructure. As a caveat, it is worth point- 
ing out that unless the cusp of the Shapiro delay is clearly 
visible in the residuals, the fit is likely to become degener- 
ate in the presence of noise, as fast-transiting small masses 
with large impact parameters will be difficult to discern from 
slow-moving massive objects with small impact parameters. 



4 DARK MATTER DISCRIMINATION 

The small-scale structure produced in our universe is very 
sensitive to the type (thermal, non-thermal, or baryonic) 
and particle properties (masses and couplings) of dark mat- 
ter. While baryons cannot be all of the dark matter, they 
do compose a significant amount of the nonluminous mat- 
ter in the universe (~ 13 per cent). Baryons do not exist 
in isolation in small mass clumps, but rather as a part of 
very large mass structures, due to the fact that their col- 
lapse is suppressed on all scales which enter the horizon prior 
to baryon-photon decoupling (at recombination). Therefore, 
even though baryons exist in the form of diffuse gas clouds 
and compact halo objects today (in addition to stars), these 
baryons were never isolated from the non-baryonic dark mat- 
ter halos in which they are found. The signature that bary- 
onic dark matter (MACHOs) will leave in the pulsar timing 
measurements is that of a point mass source, as given by 
changes over time in equation (|8]). For negligible changes in 
I and a;||, this would result in a pulse-arrival-time shift of 

At,-AU ^ 2.95x10-" s ( M^] In ( + -] 

between n pulses. Additionally, MACHOs are compact 
enough that effects such as gravitational microlensing ought 
to be observable if one knows where to look as well, allowing 
for possible cross-correlation effects. 

The vast majority (~ 87 per cent) of dark matter is 
non-baryonic, however (contained in Odm). In contrast to 
baryons, the non-baryonic matter is practically collision- 
less, and thus collapses to form much more diffuse struc- 
tures than ordinary matter. Additionally, the suppression 
of non-baryonic structure ceases at a much earlier epoch, 
as thermal dark matter kinetically decouples from the pri- 
mordial plasma at a typical temperature of 0(10 MeV), and 
non-thermal dark matter is always decoupled both thermally 
and kinetically from the plasma. The exact epoch of kinetic 
decoupling of thermally produced dark matter is determined 
by the mass of the decoupled particles and its interactions; 
this in turn dete rmines the mass fun ction of small-scale dark 
matter clumps (jBertschingerl [200^ . As a result, thermally 
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produced dark matter will have a sharp drop in its mass 
function for masses below ~ 20 A/©, with the exact value 
dependent upon the exact particle properties chosen. Non- 
thermal dark matter, by contrast, will not only have no such 
feature, but may also have a sharp enhancement in the num- 
ber density of collapsed structures on very small scales under 
the right conditions. For the case of non-thermally produced 
axions where the Peccei-Quinn symmetry is broken subse- 
quent to the end of inflation, the initial density fluctuations 
on scales correspond ing to the h orizon at T ~ Aqcd will 
be greatly enhanced llSikividbOOel ). This will result in much 
more numerous and denser structures on these small scales 
(M ~ Mq) as compared with thermal dark matter, al- 
though they are still quite diffuse as compared to baryons. 
While mass scales of 10"'' M© are far beyond the reach of 
any foreseeable technology using pulsar timing, this density 
enhancement will affect larger mass scales as well, and may 
be detectable if there is a substantial enhancement at mass 
scales which can be probed, either today or in the future, 
by pulsar timing. 

The surefire way to distinguish between the various 
types of dark matter is to probe the low end of the mass 
function, where WIMPs and NACHOs differ from one an- 
other in their abundances. Thermal relic dark matter begins 
to experience a sharp dropoff in the number of collapsed 
structures at masses below ~ 20 M© , although that number 
is model dependent. If the abundances of these dark mat- 
ter microhalos can be measured with sufficient accuracy at 
small mass scales, it may be possible to reconstruct the low 
end of the mass function, thus constraining the types and 
properties of dark matter. It is also of interest to probe the 
density profiles of these structures, which may shed light 
on not only whether axion miniclusters are present in the 
universe, but also on the nonlinear collapse of small-scale 
structure. 

In order to calculate the types of observations neces- 
sary to probe the effect of transiting dark matter, we must 
calculate the expected transit rate of these dark microha- 
los. Moreover, as it is the mass function we desire to probe, 
different halo masses m ust be considered and sea rched for. 
Press-Schechter theory (jPress fc Schechterl Il974l ) indicates 
that the number density of collapsed structures which should 
be found in our universe today is strongly dependent upon 
the mass of the structure of interest. At small scales, the 
mass function is more accurate if replac ed by the Sheth- 
Tormen function (jSheth fc Tormenlll999l ). The mass func- 
tion dn/dM then obeys the equation 



M dn _ d(\nv) 
d In M ^ d (In M) 

with 



vf{v) = 0.322 



1 



V2 



V2v 



2V2 



(13) 



(14) 



where v = cr /5c, 5c — 1.686, and a is the root-mean-squared 
mass density perturbation in a given region of space. The 
mass density perturbation spectrum is highly dependent 
upon the transfer function and power spectrum of the dark 
matter, and to a lesser extent, t h e win dow function. It is 
estimated in iLoeb fc Zaldarriagal (|2005l ) that there should 
be a significant fraction of the Milky Way halo in dark mat- 



ter clumps of approximately 10'^ Af©, implying that there are 
10" to 10^^ of these clumps in our galaxy. (A more accurate 
estimate of the mass function, although it assumes a specific 
model of thermal dark matter, can be found in iBertschingeil 
(j2006t ).) Assuming a pulsar-to-observer distance of 2kpc, a 
typical transverse velocity of a dark matter microhalo of 
vi_ ~ 200 km s~^, 10^^ microhalos of mass lO'^ A^©, and an 
NFW profile for the dark matter structure in our own galaxy 
with a turnover radius of 25 kpc, the local number density 
of dark matter microhalos is 



PmH 



10-4 M, 







133 pc 



(15) 



based on a p(r) for the Milky Way given by equation (|ll|l 
with Tmax = ^bro and our radius from the center of the 
galaxy taken to be 8 kpc. We note that this estimate is at 
the high end of the possible number density of microhalos, 
differing by as much as a factor of O(IO^) from other esti- 
mates. This results in a transit rate near the line-of-sight to 
a pulsar of 



= 0.034 



133 pc" 



V200 kms-iy 



10-* pc 
transits 



I 



year 



2 kpc 
(16) 



where bmax is the maximum impact parameter a microhalo 
can have relative to the line-of-sight and still be of interest 
on our observable (~ 25 years) timescales. This rate cor- 
responds to approximately one microhalo with the above 
properties transiting over the observed timescale of the pul- 
sar. The transit rates for different mass scales and number 
densities can be computed in a similar fashion for differing 
mass functions. 

The estimates given here are for dark matter microha- 
los embedded in the galactic halo. The presence of collapsed 
baryonic substructure, however, particularly in the galactic 
disk, may alter this profoundly. Extended, low density as- 
trophysical objects such as dark matter microhalos may be 
gravitationally disrupted (or even captured) by much more 
dense, massive baryonic st ructure in the disk as the low den- 
sity objects pass through (|Goerdt et al.|[2006l ). On the other 
hand, this may not be the globular clusters, which 

are somewhat more extended and significantly more mas- 
sive, can remain in tact and appear to be undisturbed by the 
galactic disk (.Kobulnickv et al]|2005l ). It is thus unknown 
whether dark matter substructure will survive intact in the 
galactic disk. There may therefore be a tremendous amount 
of information obtainable by comparing the time-delay ef- 
fects for pulsars found in the galactic halo, where dark mat- 
ter substructure may be relatively undisturbed, with those 
found in the galactic plane. 

The net signal which ought to be searched for in the pul- 
sar data is a superposition of all the baryonic MACHOs and 
dark matter microhalos which transit a significant distance 
in the x± direction with small impact parameters relative to 
the 1.0. s. of the pulsar of interest. Based on the transit rate 
for low-mass dark matter microhalos computed in equation 
(|16|) . it ought to be a rare event to observe a transit with a 
large enough mass (~ 10* M(s) to be detectable. With 0(10) 
presently observed pulsars which are sufficiently stable over 
their observed timescales of 0(10 years), there is a small but 
reasonable probability that evidence for dark matter micro- 
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halos exists in the current pulsar data. Although there are 
over 100 millisecond pulsars presently known, many of them 
are either weak radio sources or located in globular clusters, 
rendering them incapable of detecting dark matter through 
precision timing campaigns with current instrumentation. 
The transit rates are small enough so that spatial averaging 
is not necessary, as only a small number of these dark matter 
MACHOs (which behave as point sources) and microhalos 
(which behave as extended sources) of sufficient mass should 
be observable at all. 



5 DISCUSSION 

This paper has demonstrated the feasibility of using millisec- 
ond pulsars to detect dark matter within our own galaxy. 
Although the dark matter is distributed in a smooth halo 
around our galaxy to a zeroth approximation, departures 
from a smooth density distribution (i.e., dumpiness, cuspy- 
ness, and caustics) are features of realistic models of dark 
matter. The dark matter substructure found within our 
galaxy can be quantified through theoretical modelling of 
small-scale structure through the deeply nonlinear regime, 
but can also be probed experimentally through precision 
measurements of pulsar timing. As the dark matter distri- 
bution along the line-of-sight to a pulsar changes over time, 
the light-travel-time changes as well (as given by equation 
(0} for an individual source). This effect can be quantified 
for each source, its effect on the pulse arrival time can be 
calculated as detailed in Section |3l and summed over all rel- 
evant transiting sources. We have illustrated that the rate 
of interesting events is non-neligible, and such a signal may 
already be present in archival data. A successful detection 
of dark matter substructure using pulsar timing would be 
the first definitive detection of dark matter within our own 
galaxy. 

Beyond merely detecting dark matter substructure, pul- 
sar timing has the potential to provide information about 
the masses and density profiles of dark matter substruc- 
ture. Probing galactic substructure is the only method at 
present capable of measuring the mass function on scales 
below ~ 10® Mq . Accurate measurements of the mass func- 
tion on small scales can provide insights into dark matter 
particle properties such as the type of dark matter (WIMPs 
or NACHOs), the temperature of kinetic decoupling, and the 
dark matter's mass and scattering cross sections. Addition- 
ally, accurate measurements of the density profile of these 
dark matter substructures, which is directly related to the 
time delay signal to be probed by pulsar timing, provides 
insights into both the particle properties and the nonlinear 
evolution of dark matter. The density profiles may be signif- 
icantly different than standard cold dark matter simulations 
indicate due to dark matter substructures interacting with 
one another and with baryons, and due to the pr ocess of 
tidal stripping, as analyzed in (jCoerdt et al ] l2006h . While 
the survivability of these structures is still under debate, it 
appears to be more of a quantitative question (i.e., what are 
the percent of surviving halos, what percent are disturbed 
and to what level, how much of the denser cores survive) 
than a qualitative one. If strong limits can be placed on the 
presence of dark matter substructure in the galactic plane 



from pulsar timing, it may be evidence indicating that tidal 
disruption of these small-scale structures is very efficient. 

In addition to the effects due to intervening dark mat- 
ter, there are two other interesting physical phenomena 
which can i nterfere wit h this signal. A gravitational wave 
background I Jenet et al.l (|2006| ) is expected to have individ- 
ual waves containing important effects at the time delay 
level of 0(1 ns), with the overall stochastic background con- 
tributing a cumulative effect of 0(100 ns). Although the sig- 
nal has specific templates which are quite different from 
that of transiting dark matter, the pulse delay is of the 
same magnitude for 0(10'^ M^) substructure as it is for a 
gravitational wave background, thus care must be taken to 
distinguish these signals. Additionally, there are baryons in 
the interstellar medium that can alter the pulse arrival time 
(|Weisberg fc St animirovic''2007^ , although their effects may 
be correctable (You ct al. 2007). Although these signals are 
well-studied, they are anticipated to have significantly larger 
effects on the pulse arrival time, and thus they must be un- 
derstood with great care in order to be properly subtracted 
out. 

The possibility of detecting dark matter using pulsar 
timing is exciting and quite realistic. With the transit rate 
given in equation (|16p and the fact that there are ~ 10 
pulsars with the necessary properties for detection, the lim- 
iting factor is only to understand the pulsars' timing resid- 
uals. Recent advances have allowed pulse-to-pulse residuals 
for some pulsars to be constrained to 0(0.1 /is) with ~ hour- 
long integration times (|van Straten et al ]|200ll : |Jacobv" et al.l 
I2OOS), corresponding to an improved dark matter mass de- 
tectability limit of ^ 10^ for these pulsars. There is an 
exciting possibility ([Edwards. Hobbs, fc Manchesteij |2006| ) 
that pulsar timing may soon be accurate to the O(lns) 
level, as the software is no longer a barrier to such accu- 
racy. Significant improvements in measurement sensitivity 
would need to be made as well, of course, but if that hap- 
pens, it would allow for the detectability of ~ 10 ob- 
jects, probing past the expected cutoff in the mass func- 
tion for WIMP dark matter. A lthough it is estimated in 
lEdwardsT Hobbs. fc Manchested ((2006) that a factor of 10 
improvement is needed in sensitivity to obtain timing accu- 
racies of ~ 1 ns, the upcoming Square Kilometer Array will 
have the sensitivity necessary to improve the residuals by 
up to two orders of magnitude for many tens of pulsars. It 
is not, therefore, overly optimistic to consider the method of 
pulsar timing as having the potential to detect dark matter. 
We also point out that the use of averaging techniques may 
allow for a reduction in random errors in the pulsar resid- 
uals. When the Square Kilometer Array is complete, and 
routinely timing hundreds of pulsars to sub-/is accuracy, it 
is possible that the presence of dark matter clumps may 
limit the achievable accuracy of pulsar timing, dependent 
upon the masses and number densities of the dark matter 
clumps. 

An accurate mass function can be derived from 
Press-Schechter and/or Sheth-Tormen theory via equations 
(|13ll4p and the linear power spectrum and transfer function, 
which can be obtained for any given model of dark mat- 
ter. Realistically, the mass function can be probed all the 
way down to this detectability limit and perhaps further, 
if superior analysis techniques are developed. We empha- 
size that probing dark matter substructure within our own 
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galaxy is the only handle on structure formation at scales 
below ~ 100 kpc at present. Furthermore, there may also be 
a correlation between dark matter microhalos that are either 
sufBciently massive , dense, or nearby a nd a possible WIMP 
annihilation sie nal (Ide Boer et al.ll2005l ). which may provide 
further information about the nature of dark matter. 

Uncovering the nature of dark matter remains one of the 
most important problems in cosmology today. Precision pul- 
sar timing measurements can be used to probe dark matter 
substructure transiting near the line-of-sight within our own 
galaxy. Changes in the light-travel-time of emitted pulses, 
although small (of order 10~® s) and over long (~ 10 years) 
timescales, are nonetheless observable in practice based on 
the current accuracy of pulsar timing. Improvements in the 
understanding of pulsars themselves, as well as the local en- 
vironments of the pulsar and our solar system, will allow us 
to probe smaller and smaller mass scales. Additionally, in- 
formation gained about the transiting substructure can be 
used to constrain both the particle properties of dark mat- 
ter and the density profiles of that structure. Overall, pulsar 
timing measurements may provide a new window into un- 
covering the properties, method of creation, and evolution 
of the non-luminous matter in our universe. 
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